arXiv:1507.02034vl [physics.flu-dyn] 8 Jul 2015 


Abrupt growth of large aggregates by correlated coalescences in tnrbulent flow 

Jeremie Bee/ Samriddhi Sankar Ray/ Ewe Wei Saw/’^ and Holger Homann^ 

^Laboratoire J.-L. Lagrange, Universite Cote d’Azur, OCA, CNRS, Bd. de I’Observatoire, 06300 Nice, France 
^International Centre for Theoretical Sciences, Tata Institute of Fundamental Research, Bangalore 560012, India 
^Laboratoire SPHYNX, SPEC, CEA Saclay, CNRS, 91191 Gif-sur-Yvette, France 

Smoluchowski’s coagulation kinetics is here shown to fail when the coalescing species are dilute 
and transported by a turbulent flow. The intermittent Lagrangian motion involves correlated violent 
events that lead to an unexpected rapid occurrence of the largest particles. This new phenomena is 
here quantified in terms of the anomalous scaling of turbulent three-point motion, leading to signif¬ 
icant corrections in macroscopic processes that are critically sensitive to the early-stage emergence 
of large embryonic aggregates, as in planet formation or rain precipitation. 
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The formation of planets in circum-stellar disks [DE] as 
well as the initiation of rain in warm clouds mm in¬ 
volve the coalescence of small dilute bodies suspended 
in a highly turbulent flow. It is crucial, in both cases, 
to determine the speed at which the largest objects are 
formed. Massive planetary embryos or big raindrops de¬ 
couple from the underlying flow and accrete smaller par¬ 
ticles more efficiently mm- They are, very likely, the 
precursors for a run-away growth and possibly trigger the 
full coagulation process. Turbulent fluctuations might 
be essential in the formation of such large objects mm 
but their precise role is still far from being fully under¬ 
stood. Significant progress has been made in understand¬ 
ing the enhancement of kinetic collision kernels due to 
turbulence. It is important to recall two key mechanisms 
present in the particle dynamics: preferential concentra¬ 
tion |I 0 j . giving rise to high densities, and the sling ef¬ 
fect [ID or caustic formation [I 2 j . responsible for large 
velocity differences; both mechanisms enhance the rate at 
which particles approach each other. Precise quantitative 
models accounting for these two effects require appreci¬ 
ating the influence of turbulence [num. However, their 
origin is not directly related to turbulent fluctuations but 
rather comes from the inertia of the suspended particles 
and the resulting detachment of their trajectories from 
the flow. Their impact on collision rates can then be 
studied in simple random flow |1514I7] . 

In this Letter we show that, by its own, turbulent 
transport speeds up the growth of large objects. In the 
Lagrangian evolution of fluid elements, scaling and ge¬ 
ometry are tied up by non-trivial memory effects. These 
interdependences lead to intermittent multiscaling prop¬ 
erties of advected passive scalar fields [T8j|T9|. In the con¬ 
text of growth by coagulation, as shown in this work, they 
are responsible for a power-law tail in the distribution of 
times between successive collisions, yielding intricate cor¬ 
relations in the sequence of coalescences experienced by 
individual particles. Because of this effect, we find that 
the number of large objects grows as a power law at short 
times, with an exponent much smaller than the one ob¬ 
tained from kinetic population-balance approaches. The 


value of this exponent is expressed in terms of the anoma¬ 
lous scaling exponent /a associated to the third-order cor¬ 
relations of an advected passive scalar. 

To simplify the presentation, we focus on an initially 
mono-disperse suspension consisting of ni monomers (T) 
with mass mi. The extension to poly-disperse situations 
is straightforward. These particles evolve in a turbulent 
flow and might coalesce, summing-up their masses, when 
they collide. This dynamics leads, after sometime, to 
the formation of a broad spectrum of particle sizes. We 
denote by (T) those constituted of i monomers and thus 
with a mass i x mi. Our goal is to determine how fast the 
number ni(t) of particles (T) grows with time for i > 1. 
Simple population-balance considerations lead to 
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^*(0 = A X! ~ X! (1) 
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where the dot denotes time derivative. Qi,j{t)dt is the 
number of coalescences (f) + Q) occurring between times 
t and t -f dt. The first term in the right-hand side, the 
source, accounts for the rate at which particles (T) are 
created. The second, the sink, handles the coalescences 
of such particles with all others. When n/O) = 0 (for 
i > 1), the global coalescence rate Qij can be written in 
terms of the individual particle rate by summing over all 
the creations of @’s 

= [ Xij{t - s\s)nj{t)hi{s)ds. (2) 

Jo 

Aij(r|s) is the rate at which a particle @, created at 
time s, coalesce with a (J) at time s -I- r. For statistically 
steady particle dynamics, this quantity is independent of 
the creation time s and Aij(T|s) = Xi^ir). Also, this 
rate relates to the probability distribution Pi,j{T) of the 
time to next collision, which is given by 

A.(T) = A,,,(r)e-/o^"m(E)dr'_ ( 3 ) 

This is the distribution of waiting time associated to the 
non-homogeneous Poisson process with rate Aij(r). 
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At sufficiently long times t, the coalescence rate 
is expected to approach a finite limit . Successive col¬ 
lisions of a single particle then appear to be uncorrelated. 
They define a memoryless process and Pijir) tends to the 
exponential distribution with rate parameter X°°. The 
population-balance system then reduces to 

^ i—1 oo 
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This is the celebrated Smoluchowski coagulation equa¬ 
tion [20]. The stationary rates are usually referred 
to as the collection or coalescence kernels. The work 
cited above on particle inertia was actually devoted to 
estimating their dependence upon particle sizes and the 
turbulent fluctuations of the carrier flow. The kinetic 
model Q leads to predictions concerning the short-time 
increase of the Ui’s. At the early stages of particle growth, 
the number ni of monomers remains almost constant and 
creations are dominant in the population balance. We 
thus get n 2 — Xfin\/2, so that n 2 {t) ~ n\X'^it/2. 
For the next size, we have ha ~ X ^2 ''^2 and thus 

n 3 {t) ~ nf A“i A “2 f^/4. We obtain recursively 

ni{t) ~ n\ , (5) 

where the times ti are averages of the times 1 /A^^, associ¬ 
ated to the different combinations of coalescences ( 2 )+® 
that are necessary to form a particle (T). The consistency 
of the assumptions can be checked a posteriori: The cre¬ 
ation terms in 0 are always (x and thus prevail at 
short times over the dominant destruction term ex 

The main assumption leading to Smoluchowski kinetics 
0 is a convergence of the coalescence rate to its limiting 
value A°° much faster than the evolution of n^. This is 
ensured for instance when the particles are very dense. 
For explaining the formation of large particles in a dilute 
suspension, these timescales are in general not sufficiently 
separated. The sudden appearance of sizable aggregates 
requires a brisk sequence of coalescences that are very 
likely to be correlated to each other. When, in addition, 
the coalescing species are transported by a turbulent flow, 
such correlations speed up the growth of large particles. 

A statistically steady turbulent flow involves interac¬ 
tions between eddies of various sizes, ranging from the in¬ 
tegral scale L, where kinetic energy is injected at a rate e, 
down to the dissipative scale 77 = below which 

viscous damping dominates {v denotes the kinematic vis¬ 
cosity of the fluid). The degree of turbulence grows with 
the extension of this spatial span and is measured by 
the Reynolds number Re = The intermediate 

scales between 77 and L define the inertial range through 
which energy cascades with a rate e. Dimensional argu¬ 
ments suggest that the velocity increments between two 
points separated by a distance r in the inertial range be¬ 
have as Ur Such a phenomenology, referred to 



FIG. 1. Left: Distance traveled by fluid elements in a 3D 
turbulent flow during one large-eddy turnover time. Long 
(white) and short (purple) distances, represented here as a 
function of the final position in a 2D slice, define an intri¬ 
cate landscape with fronts where particles coming from far 
apart meet together. Right-. Sketch of the event leading to 
correlated successive collisions. At the initial time (top), two 
particles i and k are located at a collision distance r ^ rj, 
while a third one j is at distance r' ^ rj far from them. A 
time T later (bottom) j has approached i at a distance R' g 
while k, having collided or not, has escaped to i? S> 77 . 

as Kolmogorov 1941, is often enough for capturing the 
most signiheant effects of turbulent fluctuations; in real¬ 
ity the scaling properties display slight deviations, due to 
intermittency, from this dimensional prediction |21L I22j . 

The breakdown of scale invariance is much more strik¬ 
ing for mixing statistics, owing to the fact that turbulence 
mingles together fluid elements in a robust manner. This 
pops up with the presence of quasi-discontinuities in the 
Lagrangian map where materials originating from dis¬ 
tinct regions of the flow are violently brought together 
(see Fig. Left). The emergence of such fronts is due 
to the inertial-range roughness of the velocity held and 
the associated non-uniqueness of fluid element trajecto¬ 
ries. Two initially separate tracers Xi{t) and X 2 {t) that 
closely approach each other become indistinguishable and 
separate afterwards following Richardson’s superdiffusion 
\xi{f) — X 2 {t)y ~ et^. Still, when interested in more 
than two fluid elements, this explosive behavior is con¬ 
strained by the underlying presence of statistical conser¬ 
vation laws induced by the spatial correlations of the ve¬ 
locity held [181 mi US] • There exists specific functions of 
the shape and size of a cloud of n tracers that on average 
do not vary with time. Such zero modes are known to 
yield anomalous scaling in the statistics of an advected 
passive scalar 6. Its structure functions behave in the in¬ 
ertial range as {{9{x-\-r) — 6{x))'^) ~ -v^rhere the 

discrepancies 5n of the exponents from their dimensional 
prediction relate to the anomalous scaling of the transi¬ 
tion probability of the distances between n tracers. As 
we will now see, the behavior of the three-point motion 
is in fact of relevance to coalescences. 

In dilute suspensions, a coalescence results from two 
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successive processes. First, the turbulent flow needs 
to bring two initially separate particles at a sufficiently 
close distance < r]. Second, these close particles need 
to actually merge, and this involves various microphys¬ 
ical mechanisms (particle inertia, hydrodynamical inter¬ 
actions, surface effects). This leads us to write the coa¬ 
lescence rate as a product of two contributions: 

( 6 ) 

The contribution from turbulent transport can be written 



u,^Pz{.R,'n,T\r],r', 


0) dr'di? 


(7) 


involving the transition probability pa of the three-point 
motion. More specifically, two successive collisions occur 
if three particles (see Fig.^Right), initially separated by 
distances r = i] and an arbitrary r', come in a time r to 
distances R, arbitrary, and R' = r]. The relation Q is 
obtained by integrating over all possible initial distances 
r' of the particle (J) from the particle (?), with a weight 
oc given by a uniform three-dimensional spatial dis¬ 
tribution. The three tracers Q) and undergo in 
a time r an evolution from a degenerate triangle with 
r <C r' to another degenerate triangle with, this time, 
R ^ R'. In turbulence, the probability transition be¬ 
tween such configurations can be written as 


P 3 (i?,i?',r|r, /, 




where ^3 is the anomalous part of the scaling exponent 
associated to the third-order statistics of an advected pas¬ 
sive scalar; its value is universal (independent of the in¬ 
jection mechanism) and « 0.18, as reported from several 
experimental and numerical studies |24j . In the expres¬ 
sion above the first factor comes from integration over 
angles and can be seen as a small solid-angle contribu¬ 
tion. The second factor originates from intermittency 
and gives a dependence upon the integral scale L. Phys¬ 
ically, it means that when ^3 > 0 , the closer is the third 
particle, the more likely it is to approach one of the other 
two. The last terms involve a dimensionless function 'k 
that imposes Richardson’s scaling for backward and for¬ 
ward pair evolution. This specific form of the three-point 
transition probability leads to 
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where tl = s is the large-eddy turnover time. 

Turbulent transport thus leads to a power-law depen¬ 
dence in time of the coalescence rate. Plugging this be¬ 
havior in the global coalescence rate ([^ and by using the 


population-balance equations Q, one obtains a short- 
time behavior of the number of particles (?) that reads 


n,{t) (9) 

Here the characteristic times ti are (x iTj^, with a propor¬ 
tionality constant that involves the various microphysi¬ 
cal rates of the coalescences leading to 0. For 

(53 > 0 the algebraic exponent appearing in Q is smaller 
than that obtained in (i from Smoluchowski’s kinetics. 
The intermittency of turbulence mixing is thus enhanc¬ 
ing the short-time growth by coalescence. In addition, 
the larger is the aggregate size considered, the stronger 
is this enhancement. Indeed, when i is large and t 
the formation of a particle @ requires a large number of 
correlated coalescences separated by inertial-range times 
and the population dynamics is dominated by 

In order to corroborate our theoretical predictions on 
the enhancement of coalescences by turbulent mixing, we 
have performed direct numerical simulations for the evo¬ 
lution of a dilute population suspended in a turbulent 
flow. We start from one billion inertial point-particles 
whose dynamics is given by a viscous Stokes drag: 

Xn = -—[Xn-u{x„,t)], ( 10 ) 

'^n 

where u designates the fluid velocity held. It is ob¬ 
tained numerically by a pseudo-spectral integration of 
the incompressible Navier-Stokes equation using 2048^ 
gridpoints. A large-scale forcing is applied in order to 
maintain the how in a developed turbulent state with 
i?e« 50 000. 

The particles follow the how with a time lag given by 
their individual response times r„. Each particle has a 
virtual radius a„(t) and t„ oc a^. We start from a mono- 
disperse suspension with monomers having initially all 
the same radius 0 ^( 0 ) = oq « r?/10. When two particles 
approach at a distance equal to the sum of their radii 
(detected using a billiard algorithm), they merge, con¬ 
serving mass and momentum. Particles inertia is mea¬ 
sured by their Stokes numbers S't„ = t„ which 

is initially small Stn{0) « 0.1. Inertia effects can thus be 
clearly neglected when interested in inertial-range length 
or time scales. The suspension is dilute: their volume 
fraction is approximately 5 x 10“®, which represents in 
our how one particle for each cube of volume 10 rf and is 
consistent with, for example, typical settings in a warm 
cloud of our atmosphere. 

Figure shows on log-log scales the time evolution 
of the number ni{t) of particles made from the merger 
of i monomers (that is with radius a(t) = for 

? = 1,2,... 15. Data (dots) is approximated very well by 
ni{t) oc f0.73(2-2)-i-i ((Jagiied lines), corresponding to the 
predicted power laws (|^ with (53 = 0.18. Such a behav¬ 
ior persists for times larger than the large-eddy turnover 
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FIG. 2. Time evolution of the number riiit) of particles (nor¬ 
malized by the initial number of monomers ni(0)); the mass 
i increases from 1 to 15 from top to bottom. The dots are the 
results of direct numerical simulations and the dashed lines 
show for i > 2 behaviors oc t0.73{i-2)-i-i from for 

63 = 0.18. 

time. The result of our simulation confirms the accu¬ 
racy and the relevance of the predictions made earlier in 
this Letter. Large aggregates are appearing faster than 
predicted from kinetic models. 

To confirm that this enhanced growth is indeed result¬ 
ing from correlated successive collisions, we have mea¬ 
sured the probability density pi^i(T) of the time lag r be¬ 
tween the creation of a particle (T) and its next collision 
with a (T). Results are shown in Fig. for i = 2,3,4, 5. 
The distributions clearly display for t in the inertial 
range a power-law decay, followed by a (stretched) ex¬ 
ponential cutoff at T ^ T/,. The measured value of the 
algebraic exponent is consistent with the predicted value 
— (3/2)(53 « —0.27. This confirms that the inter-collision 
time distribution follows (|p) with a time-dependent coa¬ 
lescence rate Ai.j(r) oc 



Time r to next collision (in units of re) 


FIG. 3. Probability density function of the time between suc¬ 
cessive coalescences of a particle i with a monomer of mass 
1. The power-law tail has, as we predict theoretically, an 
exponent « —0.27 —(3/2)53 


In conclusion, let us stress again that intermittency of 
turbulent mixing is responsible for an enhanced growth 
of dilute coalescing aggregates. To our knowledge, this 
is one of the first instances where turbulent anomalous 
scaling laws play a critical role in describing to leading 
order a process with practical implications. Here, only 
third-order statistics are relevant since successive binary 
collisions involves the evolution of triplets of tracers that 
form degenerate triangles. Higher-order statistics enter 
other configurations (when, for instance, two particles are 
simultaneously formed and then merge) but they give 
subleading contributions. Finally, it is worth mention¬ 
ing that the effect unveiled here might be accounted for 
by modifying kinetic models. When coarse-graining the 
population dynamics on sufficiently large timescales, cor¬ 
related successive collisions will then appear as simulta¬ 
neous multiple collisions. 
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